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Abstract 

Within perturbative QCD approach based on kx factorization, we analyze the scalar mesons 
/o(980) and /o(1500) productions in B decays. By identifying /q(980) as the composition of ss 
and nn = (uu + dd)/^/2, we calculate the exclusive decays B — > /o(980)i^. We find that the 
non-factorization /o-emission diagrams can give larger contribution to the branching ratio, than 
the previous PQCD calculation. Our new results can explain the current experimental data 
well. Under the assumption of quarkonium dominance, we study the branching ratio of decays 
B — > /o(1500)i4T. The results show that in the two-quark picture of /q meson the contribution 
from ss component is at the similar size as that from the nn component. Comparing with the 
data, our results show the preference of /o(1500) as a member of the ground state of scalar qq 
nonet. Similar results can also apply to /o(1370) and /o(1710), if these mesons are dominated 
by the quarkonium content. With more experimental data in future, these studies will help us 
understand the intrinsic characters of these scalar mesons. 



1 Introduction 

In spite of the striking success of QCD theory for strong interaction, the underlying structure of the 
light scalar mesons is still under controversy theoretically [H [2] • In the literature, there are many 
proposals such as qq, qqqq, meson-meson bound states or even supplemented with a scalar glueball. 
It is very likely that they are not made of one simple component but are the superpositions of these 
* Mailing address. 
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contents and it is model dependent to determine the dominant component. The different scenarios 
may give very different predictions on the production and decay of the scalar mesons which can be 
tested by the related experiments. Although intensive study has been given to the decay property 
of the scalar mesons, the production of these mesons can provide a different unique insight to the 
mysterious structure of these mesons, especially their production in B decays. Compared with D 
meson decays, the role of scalar particles in B decays is much more noticeable because of the larger 
phase space. 

/o(980) is the first scalar meson observed in B decays with the decay mode B — > fo(980)K. 
In the three-body decays B^ — > K^tt^ 1 ^ [3], Belle found a large branching ratio for B~ — > 
X~/o(980) — * K~ (7r + 7r~), which was confirmed by BaBar [3] later. Using the branching fraction 
of /o — ► 7r + 7r~, we can obtain a large branching ratio at order 10~ 5 for the decay B — > /o(980)i^. 
These measurements of the decay B — > fo(980)K has arisen much interest on theoretical side. The 
earlier Perturbative QCD (PQCD) approach calculation [5j [6] shows a smaller branching ratio than 
the experimental data. Recently, these modes have also been studied to probe the new physics 
beyond standard model in [9] using the generalized factorization approach. They find that in 
standard model, the branching ratio is quite below the experimental values and therefore, these 
modes can be a probe of the R-parity violation supersymmetric model. Within the framework of 
QCD factorization approach (QCDF), B — * fo(980)K has also been studied recently [TJ [8]. With 
the decay constants and light-cone distribution amplitudes derived from the QCD sum rules, the 
updated results in QCDF [8] suffice to explain the experimental data. It is necessary to re-analyze of 
these decay channels in PQCD in order to find out whether the differences arise from the difference 
between the two approaches or only the different non-perturbative inputs. 

For B — > /o(1500)i^, there is a puzzle in experiments: both Belle [TO] and BaBar [11] found 
a resonance in the K + K~ mass spectrum of B — > (K + K~)K decays, whose mass and width is 
consistent with / (1500). Due to the large ratio r(/ (1500) -> vr7r)/r(/ (1500) -> KK) = 4.06 [TO], 
we expect the similar peak in the corresponding tttt channel. But there is no signal in the decays 
of B —* K( / k + it~) |1Q|. I13j. In order to make it clear, more experimental data are required. On the 
other side we should also know the theoretical predictions onB-> /o(1500)i^. 

Minkowski and Ochs [TO] have studied B — ► /o(1500)if by the assumption that /o(1500) is 
pictured as the superposition of ss and fin ^7\. In their study, the decay amplitudes are dominated 
by the QCD penguin operators b — ► sqq (q=u,d,s) and the chromomagnetic penguin operator 0$ g , 
while the tree operators are neglected for the suppression of CKM matrix elements and annihilation 
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topology contribution is also omitted due to the power suppression of 1/mg. But it is shown that 
the annihilation diagrams are not negligible in B — > ttK etc. [18] which give a large contribution 
to the imaginary paiiQ. This implies the annihilation contribution may not be negligible in B — > 
/o(1500)X either. In this paper we perform the PQCD study on this decay mode to provide a 
systematic and reliable analysis. 

The outline of this paper is as follows: In SectJJl we briefly discuss the status of the study 
on the physical properties of /o(980) and /o(1500). In Sectj3j we calculate the decays in PQCD 
approach with discussions. The final part contains our conclusions. 

2 Physical properties of / (980) and / (1500) 
2.1 Quark Structure 

Although the quark model and QCD have achieved great successes, the inner structure of scalar 
mesons is not well established theoretically. There are many scenarios for the classification of the 
scalar mesons. In the scheme proposed in [17], the a and ac(900) are not considered as physical 
states. The lowest scalar qq nonet is rather formed by the iso- vector ao(980), the iso-scalars /o(980), 
/o(1500) and the iso-doublet Xq(1430). In the second scheme, it has been suggested that the light 
scalars below or near 1 GeV- / (600) (or a), / (980), i^(800) (or k) and a (980)-form an SU(3) 
flavor nonet, either qq or qqqq, while scalar mesons above 1 GeV, namely, /o(1370), ao(1450), 
Xq(1430) and /o(1500)//o(1710), form another nonet. According to the different descriptions for 
the first nonet, this scheme is divided into two different scenarios further, which we will denote as 
scenario I and scenario II respectively in this work. 

In scenario I, the first nonet is viewed as qq states. In this scenario, /o(980) is mainly an ss 
state and this is supported by the data of Df — > foir + and (ft /o7- However, there also exists 
some experimental evidences indicating that /q(980) is not purely an ss state. The observation of 
T(J/i/j — > foul) ~ nT(J/ift — > fo<ft) [E] indicates the existence of the non-strange and strange quark 
contents in /o(980). Therefore, /o should be a mixture of ss and fin = (uu + dd)/y/2: 

|/ (980)) = \ss) cos 6 + | ran) sin0, (1) 

with 8 is the mixing angle. Experimental implications for the mixing angle have been discussed 

in detail in Ref. [20J: 9 lies in the ranges of 25° < 9 < 40° and 140° < 9 < 165°. In scenario II, 

1 In the recent study on the factorization property for the annihilation contribution using the effective theory |19] . 
the leading contributions of order a s (mj)AQCD/mg are factorizable and real. 
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/o(980) is described by a four-quark state, which is too complicated to be studied in a factorization 
approach. In order to give quantitative predictions, we work in the scenario I for /q(980) only and 
identifying it as the mixture of ss and fin. 

In both scenario I and scenario II, /o(1500) could be treated as a qq state, either the first-excited 
state or the ground state. But the case becomes complicated by the possible existence of glueball 
content. Glueball is the prediction of QCD, but any explicit evidence for a pure glueball state has 
never been confirmed in the spectroscopy of isoscalar mesons. Lattice QCD studies [23] suggest the 
mass of lightest scalar glueball lies at 1.5 ~ 1.7GeV. Among the established resonances with the 
quantum numbers to be a scalar glueball, the three mesons, /o(1370), /o(1500) and /o(1710), are 
the most natural candidates. Actually, it is likely that they are the mixtures of qq and glueball. 
Different mixing mechanisms for these mesons were proposed in the literature [17} 124"! 1251 [26| 127] . 
In different mixing mechanisms, /o(1500) is described differently which will affect the production 
in B decays. In the future, the B decay experimental data can help us to specify the right mixing. 
In the following, we assume the /o(1500) meson is dominated by the quarkonium content, i.e. 
|/o(1500)) = cos#|ss) + s'm9\nn) and leave the contribution from glueball content for future study. 

2.2 Decay constants and Light-Cone Distribution Amplitudes 

In two-quark picture, the decay constants for scalar meson S are defined by: 

{Sip) 1^27^1 10> = fsPp, (S(p)\q2qi\0)=m s fs. (2) 

Due to the charge conjugation invariance, the neutral scalar mesons /o(980) and /o(1500) cannot 
be produced by the vector current. Thus fs = 0. Taking the mixing into account, the above 
definition is changed to: 

</ n M0) = ^m f Jl, (f£\ss\0) = m f J%. (3) 

Using the QCD sum rules method, the decay constants /? and /i of /o(980) have been determined 
separately but with no great difference [7]. So the assumption of = /! works well. In the 
following, we will denote them as ff . 

The twist-2 and twist-3 light-cone distribution amplitudes (LCDAs) for different components 
of /o are defined by: 

+m f M + 4--^ S)T (x)} J u (4) 
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where n+ and ra_ are light-like vectors: n+ = (1,0, 0x),w_ = (0, 1,0t). The normalization can be 
related to the decay constants: 

1 dx^' s) (x) = f 1 dx$>f s)T (x) = 0, f 1 dx^' s)S (x) = — (5) 



In the following, we assume the SU(3) symmetry and use the same LCDAs for ss and nn. The 
twist- 2 LCDA can be expanded in the Gegenbauer polynomials: 

= ^ 7 =//(M)6 a ;(l-x)^i3 m ( / i)C^ 2 (2x-l). (6) 
» c m=l 

The decay constants and the Gegenbauer moments B m (n) for twist-2 distribution amplitude have 
been studied in [TJ E] using the QCD sum rules approach: 

Scenario I: / /o(980) = (0.37 ± 0.02) GeV, / /o(1500) = -(0.255 ± 0.03) GeV, 

^(980) = -0.78 ± 0.08, £ 3 (980) = 0.02 ± 0.07, 

Bi(1500) = 0.80 ± 0.40, 5 3 (1500) = -1.32 ± 0.14, (7) 

Scenario II: / /o (i 50 o) = (0.49 ± 0.05) GeV, 

5i(1500) = -0.48 ± 0.11, 5 3 (1500) = -0.37 ± 0.20, 

where the values for Gegenbauer moments are taken at scale /x = 1 GeV. As for the explicit form 
of the Gegenbauer moments for the twist-3 distribution amplitudes $^ and <I>J, they have not been 
studied in the literature, so we adopt the asymptotic form: 

= 27Wj>- *' = 27Wj' (1 - 2 ^ (8) 

In the previous PQCD study on B — > /o(980).?C [5], the author neglected the twist-2 contribution 
but only used the twist-3 distribution amplitude &'J(x) and proposed the following form: 

= ^={3(1 - 2xf + £(1 - 2x) 2 [C 2 3/2 (l - 2x) - 3] + 1.8C 4 1/2 (1 - 2*)}, (9) 

with C\ /2 {y) = (35y 4 - 30y 2 + 3)/8, Cl /2 {y) = 3/2(5y 2 - 1). The decay constant f f = 0.2 GeV 
was used which is close to the earlier QCD sum rules study: ff = 0.18 ± 0.015 GeV [29]. The 
parameter £ was chosen as: £ = 0.3 — 0.5. While in Ref. [6], the twist-2 distribution amplitude was 
also included: 

$ f (x) = ^=G[6x(l - x)Cl l2 {\ - 2x)], (10) 

where G ~ 1.11 obtained from the corresponding value in ao(980) given by |30j. But the decay 
constant ft = 0.2 GeV is much smaller than the recent QCD sum rules results in Eq.0. The 



reason of the difference is that the scale dependence of ff and the radiative corrections to the 
quark loops in the operator product expansion series is taken into account in Ref. [8] . The larger 
decay constant can surely enhance the branching ratio and we expect a much larger branching ratio 
for the decay B — > /o(980)iT also in PQCD approach. 

3 Calculation of decays in PQCD approach and discussions 

In the standard model, the effective weak Hamiltonian mediating flavor-changing neutral current 
transitions of the type b — > s has the form: 

10,77,89 



H*ff ~ 



-7=[ E v pb v; a (c 1 o p 1 +c 2 o^-v tb v t * £ dOi], (11) 

p=u,c i=3 



p=u,c i=3 

where the explicit form of the operator Oj and the corresponding Wilson coefficient C% can be found 
in Ref. [31]. V p ^ b , V p u) s are the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. 

In the effective Hamiltonian, the degrees of freedom heavier than b quark mass m& scale is 
included in the Wilson coefficients which can be calculated using the perturbation theory. Then 
the left task is to calculate the operators' matrix elements between the B meson state and the final 
states, which suffers large uncertainties. Nevertheless, the problem becomes tractable if factoriza- 
tion becomes applicable. The PQCD approach is one of the standard factorization approaches in 
hadronic B decay studies [18J. In this approach, we apply the kx resummation to kill the end-point 
singularities and threshold resummation to smear the double logarithmic divergence from the weak 
corrections, which results in the Sudakov form factor S and the jet function J respectively. Then 
the decay amplitude can be factorized into the convolution of the wave functions and the hard 
kernel in the following form: 

A=$ B ®H® J®S®$ Ml ®$M 2 - (12) 

The hard part H can be calculated perturbatively, while the LCDAs $e, &Mi an d <&m 2 > although 
non-perturbative in nature, are universal for all modes. They can be determined by other well 
measured decay channels to make predictions here. For example, the corresponding light-cone 
distribution amplitudes of B and K mesons are well constrained by the B — ► Kit and B — ► tttt 
decays [18]. 

The leading order Feynman diagrams are given in Fig. Q] and [2] for the ss and nn component of 
/o respectively. The decay rates of B — > JqK can be expressed as: 



Table 1: Input parameters used in the numerical calculation 


Masses 


m fn(980) = 0.98 GeV mff = 1.7 GeV, 
"»/o (1500) = L5 GeV m b = 5.28 GeV 


Decay constants 


/b = 0.19 GeV f K = 0.16 GeV 


Life Times 


t b± = 1.671 x 10- 12 s t b o = 1.536 x 10~ 12 s 


CKM 


V tb = 0.9997 V ts = -0.04, 
V us = 0.2196 Kfe = 0.00367e~ i60 ° 



in which r/ = mj/niB- A is the decay amplitude of 5° — ► /o-ff and A is the decay amplitude 
can be written as 



A ( ) = ,4^ ) x cos 6* + — x sin i 



V2' 



(14) 



with 



A: 



A 



-nn 



= -v tb v* 



F? P (a & - \a 8 ) + M L e L {C z + C 4 - - ±C 10 ) + M PR (C 5 - X -C 7 ) 



+M S e P (C 6 - l -C 8 ) + F a LL ( ai - \a w ) + F a SP (a 6 - l -a 8/ 
+M PL (C 3 - l -C 9 )+M L a R (Cs - \c 7 ) 



(15) 



A_ - 



M e LL (2C 4 + \c w ) + M! P (2C & + ±C 8 ) 



V ub V: s Me{C2) ~ V tb V* 

+F e LL '(a 4 - \a w ) + F e SP \a 6 - ^a 8 ) + M LL '(C 3 - l -C 9 ) + M LR '(C 5 - \c 7 ) 
+F a LL '(a 4 - \a w ) + F a SP '(a 6 - la*) + M L a L '(C 3 - \c 9 ) + M L a R '(C 5 - \c 7 ) 



v„. h v:, 



F LL (ai)+M^(Ci) 



LL, 



2 

v tb v* 



(16) 



if p (a 6 - ±a 8 ) + M L e L (C 3 + C A - X -C % 



-C 
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+M ip (C 5 - ^C 7 ) + Mf p (C 6 - ic 8 ) + F a iL (a 4 + a 10 ) + if p (a 6 + a 8 ) 



+A^ L (C 3 + C 9 ) + A^ P (C 5 + C 7 ) 



Li?/ 



(17) 



M L e L {C 2 ) + F e LL '(ai) + M L e L \C{) + F a LL '(ai) + M^Ci) 



M e LL (2C 4 + ^Cio) + Mf p (2C 6 + l -C 8 ) 



+F e LL '(a 4 + oio) + if P/ (a 6 + a 8 ) + M L e L ' (C z + C 9 ) + A^ P '(C 5 + C 7 ) 
+Ff L '(a 4 + a 10 ) + if P/ (a 6 + a 8 ) + M^'(C 3 + C 9 ) + A^ p, (C 5 + C 7 ) 



(18) 
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where the combinations of the Wilson coefficients are defined as usual [32] : 

ai = C 2 + Ci/3, a 3 = C 3 + C 4 /3, a 5 = C 5 + C 6 /3, a 7 = C 7 + C 8 /3, a 9 = Cg + C 10 /3, (19) 
a 2 = d + C 2 /3, a 4 = C 4 + C 3 /3, a 6 = C 6 + C 5 /3, a 8 = C 8 + C 7 /3, o 10 = do + C 9 /3. (20) 

The explicit amplitudes, for the factorizable /o-emission contribution F e (the first two diagrams 
in the first line of Fig. Q] and the first two diagrams in the third line of Fig. [2]) and nonfacorizable 
contribution A4 e (the last two diagrams in the first line of Fig. [Hand the last two diagrams in the 
third line of Fig. [2]), the factorizable annihilation F a (the first two diagrams in the second line of 
Fig. Q] an d the first two diagrams in the second line of Fig. [2]) and non-factorizable annihilation 
contribution M a (the last two diagrams in the second line of Fig. [Hand the last two diagrams in 
the second line of Fig. [2]), for the factorizable i^-emission contribution F' e (the first two diagrams 
in the first line of Fig. [2]) and nonfacorizable contribution A4 e (the last two diagrams in the first 
line of Fig. [2]), are given in the appendix [XI 

For the numerical calculation, we list the input parameters in Table [H 

3.1 Branching ratios 

At first, we give the results of the form factor Fq ~^f°^ nn ^ a t maximally recoiling: 

F B^ /0 (980) = Q47) (21) 

F B-*/o(i500) = _ a39 ScenarioIi (22) 
F s^/o(i500) = g6 ScenarioIL (23) 

These form factors are large, because the decay constants of the scalar mesons are very large. The 
minus sign of the B — > /o(1500) form factor in scenario I arises from the decay constant of /o(1500). 
If /o(980) is purely composed of ss, the branching ratios of B — > fo(980)K are: 

B(B° - / (980)K°) = (22^+2+J) x iQ-6, ( 24 ) 
B(B~ -> f (980)K-) = (24^+3+J) x 10- 6 , (25) 

where the uncertainties are from the decay constant of /o(980), the Gegenbauer moments B\ and 
B 3 . If /q(980) is purely composed of nn, the branching ratios for B — > fo(980)K are: 

B(B° / (980)K°) = (28^+6+2) x 10- 6 , (26) 
B(B~ - / (980)i<:-) = (34±1±1±§) x 1(T 6 , (27) 
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Table 2: Decay amplitudes for B° f {980)K° (xlO~ 2 GeV 3 ), where "[5]" denotes the results 
using the LCDAs proposed in [Sj, "This work" denotes the results using Gegenbauer moments of 
twist- 2 distribution amplitude from QCD sum rules and asymptotic form of twist-3 distribution 
amplitudes. 



ss 


F e 


M e 


F a 


M a 


m 


6.2 


~ 


-2.6 - 2.1? 


0.04 - 0.06i 


This work 


11 


-4.0 - 5.0* 


0.9 - 8.32 


0.31 + 0.49i 


ran 


Mj 


M e 


F a 


M a 


® 








-3.1 - 3.1? 


-0.42 + 0.13i 


This work 


52 + 55i 


—10 — 13* 


1.0 - 9.9i 


-0.38 - 0.18* 


ran 


F' 

e 


M> 






m 


-7.9 


-0.01 - 0.03i 






This work 


-2.2 


0.17 + 0.58i 







where the uncertainties are from the same quantities as above. Our results are larger than the 
earlier PQCD results [5], with the branching ratios for purely ss component 

B{B° -> / (980)if°) = 1.39 x 10~ 6 , (28) 
B(B~ -> /o(980)K~) = 1.57 x 10" 6 , (29) 

and for purely nn component, B(B° -> / (980)K°) ~ B(B~ -> / (980)iT _ ) ~ 5 x 10" 6 . 

In order to find the sources of the difference, we list the numerical results for different topology 

diagrams of B° -> / (980)K° in Table M B In the table, F e (F a ) and M e (M a ) denote as the f 

emission (annihilation) factorizable contributions and non-factorizable contributions from penguin 

operators respectively, while F' e and M' e are the contributions from the K emission diagrams in nn 

component of /o(980). denotes the /o emission non-factorizable contribution from tree operator 

O2. From the table, we can see that for ss component of /o(980) the new result is quite larger 

than the previous one [5]. The main reason is that the new decay constant for /o(980) is twice as 

the previous one. Furthermore, in Ref. [5] the non-factorizable contribution is very small, which is 

2 We reproduce the result using the previous light-cone distribution amplitudes, with slightly different convention: 
we didn't factor out the —M% for decay amplitudes; we didn't factor out the decay constant in the factorizable 
amplitudes. 
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(e) (f) (g) (h) 

Figure 1: The leading order Feynman diagrams for B° — > /o(ss)ET° decay in PQCD approach 

understandable from the amplitude for this contribution: the third diagram and the fourth diagram 
in the first line of Fig. Q] cancel with each other due to the symmetry of X2 — * 1 — X2 in the twist-3 
distribution amplitudes. But after including the twist-2 distribution amplitude for /o(980), on the 
contrary, the two diagrams do not cancel with each other due to the antisymmetry of the twist-2 
distribution amplitude and give a large contribution. Although the normalization of &f(x) 

is zero, this contribution can enhance the total branching ratio further. In the factorizable part 
of annihilation diagrams (FigfjTe) and (f)), the distribution amplitude $> s in [5] can give a large 
contribution both to the real part and imaginary part. But here we include the twist-2 distribution 
amplitude and use the asymptotic form of twist-3 distribution amplitudes. The real part becomes 
positive with a value half of the one in [5], and the imaginary part is 3 times larger. This suggests 
that the annihilation type amplitude is quite sensitive to the shape of the distribution amplitudes. 

For B — ► /o (980) (hn)K, there are two kinds of emission diagrams: /o(980)-emission (last row in 
FigE]) and i^-emission (first row in Fig|2]). In Ref. [5], both factorizable and non- factorizable con- 
tribution from the /o(980)-emission diagrams are zero. But after including the twist-2 distribution 
amplitude &f(x), the non- factorizable two diagrams do not cancel with each other, thus can give 
large contribution to branching ratio. For if-emission diagrams, although we use a larger decay 
constant for /o(980), our result is even smaller than the previous one. The reason is that after we 
include the twist-2 distribution amplitude and use the asymptotic form for twist-3 $p $T, there 
exist cancelations between the (V — A)(V — A) operators and (S — P)(S + P) operators. From the 
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(e) (f) (g) (h) 




(i) 0) (k) (1) 

Figure 2: The leading order Feynman diagrams for B° — > fo(nn)K° decay in PQCD approach, 
where q denotes u or d quark in the last row 

Table [21 we can see that the /o- em i ss i° n diagrams and the factorizable annihilation diagrams give 
the largest contribution, therefore the new branching ratio is five times as the one in Ref. [5] for 
inn) components. 

Our results are also larger than the ones in QCD factorization approach [8]: B(B~ — > f$K~) ~ 
18 x 10~ 6 for ss component and B(B~ —* /qK~) ~ 1 x 10~ 6 for ran component. The main reason 
is that in QCDF the /o-emission nonfactorizable diagrams are very small. 

In Fig. [3l we plot the branching ratios as functions of the mixing angle 8. Using the above 
mentioned range of the mixing angle, we obtain: B(B° — > foK°) = (32 ~ 36) x 10~ 6 , B(B~ — > 
f K~) = (35 ~ 39) x 1(T 6 for 25° < 9 < 40° and B{B° -> f K°) = (13 ~ 16) x 10" 6 , B{B' -> 
foK~) = (16 ~ 18) x 10~ 6 for 140° < 8 < 165°, where only the central values of other input 
parameters are used. The averaged experimental data obtained by heavy flavor averaging group 
[33] are also shown in Fig. [3j 

B(B~ ^ f (?8Q)K-) = (17.110) x 10-6, ( 30 ) 
11 



B(B° -► / (980)K°) = (11.1 ± 2.4) x 10" 



(31) 



where we find the PQCD approach results in the range of 140° < 6 < 165° suffice to explain the 
large experimental data. Thus it does not need the existence of new physics. 




80 100 120 140 160 180 

e(degree) 

Figure 3: The dependence of the branching ratios for B — > /o(980)X on the mixing angle 9 
using the inputs derived from QCD sum rules, where the doted (solid) curve is for B° — > foK° 
(B~ — > foK~). The horizontal band within the doted (solid) lines shows the experimentally allowed 
region of B° — > /o-ST (-B - — ► foK~) within 1 a error. 

Now we turn to B — > / (1500)lf decays. The branching ratios in scenario I are: 



6(5° - /o( S -s)K°) 
h{ss)K~) 
B(B° -+ f (nn)K°) 
B(B~ - fo{nn)K-) 

while in scenario II, the results are: 

6(5° - / ( S - S )K°) 

6(5- - / ( S - S )r) 

6(5° - / (nn)K°) 
- h{nn)K-) 



(4.51^) x 10- 6 , 
(4.01^) x 10- 6 , 

(5- 2 ±i:?) x io~ 6 , 

(9.2ll;g) x 10- 6 , 



^l 1 ) x 10 



6 

x 1(T 6 , 

6 



(54^ 
(56111) 

■±S°) 

(60±if) x 1(T 6 , 



(46^ u ) x 10 



(32) 
(33) 
(34) 
(35) 

(36) 
(37) 
(38) 
(39) 



where fo(ss) denotes the ss component of /o(1500) and similar for fo(nn). The uncertainty is 
coming from the decay constant of /o(1500). The decay constant in scenario II is larger than that 
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in scenario I, so we can get a larger branching ratio in scenario II for both ss and nn component, and 
the contributions from the two kinds of components are very close. In scenario I, B~ — > /o(nn)A~ 
is larger than B° — ► fo(nn)K , which ia from the large tree contribution in B~ — > fo(nn)K~ . 
This also implies that there is a large CP asymmetry in I? 1 * 1 — ► /o(nn)A . As discussed above, 
/o(1500) may be the mixture of ss and nn. So we plot the branching ratios as a function of 
mixing angle in Fig. 2] for scenario I and Fig. [5] for scenario II. For the comparable contribution 
from the nn and ss components, the variation range according to the mixing angle are not very 
large, which can be seen from FigSJand FigHJ Using the mixing mechanism for /o(1500) in [27J: 
|/o(1500)) = — 0.54|nn) + 0.84|ss) + 0.03|G) and neglecting the small component of glueball, we 
get: B(B° -> / (1500)K°) = 8.7 x 1(T 6 and B(B~ -> /o(1500)A~) = 10 x 1(T 6 in scenario I; 
B(B° -> / (1500)i?°) = 42 x 10" 6 and B(B~ -» / (1500)if _ ) = 55 x 1(T 6 in scenario II. With the 
experimental data in [10J and 8(/ (1500) -> K+K~) = 0.043, we have: 

B(B~ -» /o(1500)A^ ) = (471.9 ± 51.3) x 10" 6 , Solution I, (40) 
B(B~ -» /o(1500)A~) = (67.1 ± 14.4) x 10" 6 , Solution II, (41) 

where we identify the resonance as /o(1500). We can find that the second experimental solution 
is more appropriate which is also consistent with [16]. The experimental data of solution II for 
B~ — > /o(1500)A'~ is also shown in Fig. [Hand Fig. [5] if the resonance can be viewed as /o(1500). 
In scenario I, our central value is out of the experimental range with 3a, while in scenario II, with 
possible mixing the experimental data can be well explained. 

In the above calculation, the non-perturbative inputs for /o(1500) are the decay constants and 
the light-cone distribution amplitudes derived from the QCD sum rules for the first excited scalar 
state or the ground scalar state with the mass around 1.5 GeV. Similar results can also be applied 
to /o(1370) and /q(1710), if these mesons are dominated by the quarkonium content. 

It should be more interesting to include the contribution from glueball component in these 
decays. Thus we should consider the typical diagram in Fig. [6] and other corresponding diagrams 
due to the different emissions of the gluons, but the others are suppressed as argued in [37]. The 
complete calculation of these diagrams requires the derivation of Sudakov form factor and jet 
function for a glueball state, which is quite technical. We leave this work for future study, although 
the glue component is argued to play an important role in /q. 
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Figure 4: The dependence of the branching ratios for B — > /o(1500)-ff on the mixing angle # in 
scenario I, where the doted (solid) line is for B° — ► /o-^ (-B - — * foK~). The horizontal band 
within the solid lines shows the experimentally allowed region of B~ — > foK~ within 1 a error. 

3.2 CP asymmetries 

The results of the direct CP asymmetries are listed in Tabled In B° -> f (ss)K°, there is no tree 
contribution at the leading order, so the CP asymmetry is naturally zero. In B~ — > fo(ss)K~ , 
the tree contribution is from the annihilation diagrams which suppressed by 1/mg, thus the direct 
CP asymmetry is small. As we have discussed, the /o- em i ss i° n non-factorizable diagrams not 
only give large penguin contributions but also to the tree contributions, so the CP asymmetry 
of B~ — ► /o (980) (nn)K~ is large. The direct penguin contribution to B° — ► /o(1500)(nn)-^ in 
scenario I has some cancelation between emission diagrams and the annihilation diagrams, thus 
there is a large direct CP asymmetry. But in scenario II, the annihilation penguin diagrams 
enhance the emission diagrams, so the direct CP asymmetry in scenario II is rather small. In 
B~ — > /o(1500) {fin) K~ , the CP asymmetries are large in both scenarios. The different sign is due 
to the sign of the Gegenbauer moments. 

Now we discuss the CP violation in the neutral B° decays in which there are both direct CP 
asymmetry Af^'p and mixing-induced CP asymmetry A™p. The time dependent CP asymmetry 
of B decay into a CP eigenstate / is defined as: 

A CP (t) = A d S r P {B d - /) cos(AMt) + A%f{B d - /) sin(AMi), (42) 
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Figure 5: The dependence of the branching ratios for B — > fo(1500)K on the mixing angle 6 in 
scenario II, where the doted (solid) line is for B° — > foK° [B~ — ► JqK~). The horizontal band 
within the solid lines shows the experimentally allowed region of B~ — > foK~ within 1 a error. 

Table 3: Direct CP asymmetries (in units of %): the results in the brackets are the CP asymmetries 
for fin component, the other values are for ss component. 



Channel 


Scenario I 


Scenario II 




B° -> f (98O)K° 


0(3.0) 






B- -> / (980)if- 


-1.8(-24) 






B° -> / (1500)^° 


0(24) 


0(3.5) 




B~ -> / (1500)K- 


1.2(30) 


-1.5(-27) 





with 



,,/;,■ :r> r, l A ! 2 ~ 1 Amixm . ,• , _ 2/ '" A 



^(Bd-/) = ^p^, A^F(^-/) = IT ^, (43) 

where ?? = ±1 depends on the CP eigenvalue of /, AM is the mass difference of the two neutral B 
meson eigenstates. We can also use C, S to denote the direct and mixing-induced CP asymmetry. 
(3 is the CKM angle defined as usual [T2]. If there is no tree contribution in the amplitude, the 
direct CP asymmetry is zero and A can be related to e~ 2l/3 , so the mixing induced CP asymmetry 
is proportional to sin(2/J). For ss component of /o, even with the glue component, there is no tree 
contribution at the leading order. This channel can be used to extract the CKM angle (3. For fin 
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Figure 6: One of the leading order contributions to B — > ^(glueball)^ in PQCD, other diagrams 
obtained by attaching one or both of the two gluons to any other quark lines are suppressed. 



component, there is tree contributions. But for B° — > fo(980)K° , the tree contribution is small 
and the direct CP asymmetry is only a few percent, so this mode can serve as a possible place to 
extract f3 even when taking the mixing of /o(980) into account. In scenario II, B — > /o(1500)-?T is 
similar. In scenario I, the tree contribution to B — > fo(1500)(nn)K is relatively large and the direct 
CP asymmetry is roughly 24%. Using sin(2/3) = 0.687 [12], the mixing-induced CP asymmetries 
of B — > fo(nn)K are: 



= -0.608, 

A^{B° -» /o(1500)iT s ) = -0.633, Scenario I, 

—0.642, Scenario II. 



A%?(B" -> f (980)K s ) 
A^(B° -, / (1500)iT s ) 



(45) 
(46) 
(47) 



They are not far away from — sin(2/3) 
often used: 



-0.687. On the experimental side, the parameter AS is 



AS = Affi + tan(20). 



(48) 



As we have discussed, the direct CP asymmetry in B° — ► fo(ss)Ks decay vanishes thus the 
parameter AS is zero. For fin component, the extra tree contribution makes AS deviate from 0. If 
new physics can induce the b — > s transitions and has a different phase, it would also give a non-zero 
value to this parameter. Precise measurement of this parameter and the theoretical calculation in 
the standard model can also help us to probe new physics. Our results for these three channels in 
standard model are: 



AS(B° -> f (980)K s 



0.079, 



(49) 
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AS{B° -> / (1500)iT s ) = 0.054, Scenario I, (50) 
AS(B° / (1500)ir 5 ) = 0.045, Scenario II. (51) 

The dependence of AS on the mixing angle 6 of fo is also plotted in Fig. From this figure, we 
can see that for B° -> / (980).Ks and B° -> / (1500)iT 5 decay in scenario II, there is not much 
deviation from AS = 0. But there are large deviations from AS = for the interference of ss and 
ran component of /o(1500) in scenario I around 6 = 40°. 



0.10- 




-0.25 - 

i i i i i i i i i i i i i i i i i 
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Figure 7: The 9 dependence of the AS with the solid line for B — > fo(1500)Ks in scenario I, the 
dashed for B — ► /o(1500)i£s in scenario II and the doted line for B — ► fo(980)Kg. 



3.3 Theoretical Uncertainties 

In our calculation, one of the uncertainties is from the scalar meson decay constants. These uncer- 
tainties can give sizable effects on the branching ratio, but the CP asymmetries are less sensitive 
to these parameters. There are other uncertainty sources for both branching ratios and CP asym- 
metries: 

• The twist-3 distribution amplitudes of the scalar mesons are taken as the asymptotic form 
for lack of better results from non-perturbative methods, this may give large uncertainties. 
These distribution amplitudes needs to be studied in future work. 
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• The Gegenbauer moments B\ and B% for twist-2 LCDAs of /o(980) and /o(1500) have sizable 
uncertainties. For example, the large uncertainty of B\ in scenario I and B% in scenario II 
may change the results sizably, these parameters should be constrained in future. 

• The uncertainties, from the light pseudoscalar meson and B meson wave functions, Aqcd 
and other renormalization group parameters, 0(1/ Mb) corrections and the sub-leading com- 
ponent of B meson distribution amplitude, have been systemically studied extensively in |34j : 
the uncertainty from the factorization scale choice is within 10%; the results vary by 10 — 30% 
by changing the parameter in the wave functions. 

• The sub-leading order contributions in PQCD approach have also been neglected in the 
calculation, which were calculated in Refs. [35] for B — > tttt, ttK, etc. These corrections can 
change the penguin dominated processes, for example, the quark loops and magnetic-penguin 
correction decrease the branching ratio of B —* ttK by about 20%. We expect the similar size 
of uncertainty in B — > foK decays, since they are also dominated by the penguin operators. 

• Besides, the decay amplitude suffer other power corrections which are non-perturbative in 
nature: the long distance re-scattering effect. This effect could be phenomenologically in- 
cluded in the final-state interactions [36]. But we need more data to determine whether it is 
important in B — > SP decays. 

4 Conclusion 

In this work, we re-analyze the exclusive decays B — > fo(980)K in perturbative QCD approach by 
identifying / (980) as the composition of ss and ran = (uu + dd)/\/2. The B — * fo(1500)K is also 
analyzed in PQCD approach. Our main results are as follows: 

• Using the decay constants and light-cone distribution amplitude derived from QCD sum rules, 
we find that the PQCD results can also explain the large experimental data which agrees with 
results from QCDF. 

• The non-factorizable /o-emission type diagrams can give large contributions, although the 
normalization of the twist-2 distribution amplitude for /q is zero. 

• The B — ► /o(1500)if decay is studied under the assumption of quarkonium dominance in two 
scenarios. The ran = (uu + dd)/y/2 and ss can give similar contributions. In scenario II, the 
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branching ratios are large which can accommodate with the second solution of the experi- 
mental data, but in scenario I, the predicted branching ratio is smaller than the experimental 
data. 

• The calculation of B — > fo(1500)K decays can also be applied to /o(1370) and /o(1710), if 
these mesons are dominated by the quarkonium content. 

• The mixing-induced CP asymmetries are not far away from sin(2/3) = 0.687 for B° — > 
fo(980)K s ,L and B° f (1500)K s>L . Thus these channels can provide possible places to 
extract the CKM angle j3. 
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A factorization formulae 

In this appendix, we give the factorization formulae involved in the decay amplitudes. In the 
formulae, we choose the momentum fraction at the anti-quark, thus we should use 3>/(l — x), 
$j(l — x) = $f(x) and — a?) = — $J(x) for /q. But for simplicity, we will use to denote 

— x) in the formulae. It is similar for the pseudoscalar meson K. 

In each Feynman diagram, there may be three different kinds of operators: (V — A)(V — A), 
(V — A)(V + A), (S — P)(S +P) (arising from the Fierz transformation of (V — A)(V + A) operators). 
We will use the indices LL, LR and SP to characterize the different kinds of the operators. 

If /o is emitted, the factorization formulae for the emission type diagrams are: 

/*! /*oo 

Fe P { a i) = —16TrCFm%rfff / dxidx 3 / h\db\ b^dbs <&b(xi, &i) 

Jo Jo 

x I [$£(z 3 ) + r K x 3 ($£(2:3) - &K(?aj) + 2r K <S>%(x 2 ) 

xE ei (t)h e (xi,x 3 ,bi,b 3 ) + 2r K ^(x 3 )E ei (t')h e (x 3 ,xi,b 3 ,b 1 )^ , (52) 

for the factorizable diagrams, i.e. Fig. [T](a),(b) and Fig. [2] (i),(j), and 

Me L (ai) = -32irC F m%/y / 2N c / d Xl dx 2 dx 3 b 1 db 1 b 2 db 2 $B(x 1 ,b 1 )$ f (x 2 ) 

Jo Jo 
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I [(x 2 - 1)&k(x 3 ) + r K x 3 ($%(x 3 ) - $k( x s))] E' ei (t)h n (xx,l-x 2 ,x 3 ,bxM) 

(x 2 + X 3 )<S>k{x 3 ) - TRX3 ,x 2 ,x 3 , bi, &2)}(53) 



M L e R { ai ) 




1 />00 

dx\dx 2 dx 3 / b\dh\ b 2 db 2 <&b(xi, b{] 
o Jo 



(x 2 -l)$£(x 3 ) (x 2 ) + $J(x 2 ) 



E' ei (t)h n (x 1 , 1 - x 2 ,x 3 ,h, b 2 ) x 

+r K (x 2 ~ 1) ($&(x 3 ) - *£(a*)) ($/0 2 ) + $/(x 2 ) 
r^x 3 (*&(x 3 ) + *3rC»*)) (^/(x 2 ) - */(x 2 ) 



+£'gj(t / )/in(2;i,^2,a;3,&i,^2) x 
+^2 ($£(2:3) - S/rC**)) (^/fca) - $/(x 2 ) 
+r K x 3 (*£(x 3 ) + *^(X3)) (*/(x2) + *J(X2) 



(54) 



5P/ 



32vrC F m|j/v / 2iVc jf d Xl d X 2dx 3 M&i b 2 db 2 $ B (xi, h)® f{x2)\E' ei (t) 
O2 ~x 3 - l)$£(x 3 ) + rxx 3 ($£(x 3 ) + $#(#3)) 1 - £2, x 3, &i, 62) 

Z2$£(x 3 ) ~r K x 3 ($&(x 3 ) - ^(x 3 ))l E' ei (t')h n (xt,x 2 , x 3 , b x ,b 2 ) \ , (55) 



+ 

for the nonfactorizable diagrams, i.e. Fig. Q] (c),(d) and Fig. [2] (k),(l). While if K is emitted, the 
formulae are: 



8irCFTn%fK / dx\dx2 I b\db\b2db 2 ^B{x\-,b\) 
Jo Jo 

(1 + x 2 )$/(x 2 ) - 77(1 - 2x 2 ) ($/(x 2 ) + (x 2 ))] E ei (t)h e ( Xl ,x 2 , 61, b 2 



-2rf$f(x 2 )E ei (t')h e (x 2 ,xi,b 2 , h 



(56) 



? SPr 



(a) 



rl roc 

16-7rCi?m^/ft-r^: / dx\dx2 \ b\db\ 6 2 <i& 2 $b(xi, b\) 
Jo Jo 

- <f>f(x 2 ) + r f (x 2 $J(x 2 ) - (x 2 + 2)$^(x 2 )) E ei (t)h e (x 1 , X 2,b 1 ,b 2 ) 
+2r f <Z> s f (x 2 )E ei (t')h e (x2, xi, 6 2 , 61) j , (57) 



M L e L \ ai ) 



+ 



32irC F m%/y/2Nc / d Xl dx 2 dx 3 / M&i M&2 $b(xi, fci)$£(x 3 ) 

(x 3 - l)*/(x 2 ) -r/x 2 (*/(x 2 ) - $/(x 2 ))] E' ei (t)h n ( Xl ,l - x 3 ,x 2 ,h,b 3 ) 
(x 2 + x 3 )$/(x 2 ) +r/x 2 ($f(x 2 ) + $J(x 2 ))] £^(OMxi,x 3 ,x 2 ,&i,& 3 )} (58) 
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Mz W {ai) = 32irC F m%/V2Ncr K / d Xl dx 2 dx 3 / M&i b 3 db 3 ®b{xi, &i) 

Jo Jo 

x^E' ei (t)h n ( Xl ,l -x 3 ,x 2 ,b 1 ,b 3 ) x [(x 3 - l)$/(x 2 )($£(x 3 ) + $£(x 3 )) 
+r/$J(x 2 ) ((x 2 + x 3 - 1)^(^3) + (-X2 + x 3 - 1)^(^3)) 
+r f $ s f (x 2 ) ((x 2 - x 3 + l)$£(x 3 ) - (a?2 + £3 - l)*^(x 3 )) 

- [x 3 *f(x 2 ) ($£(z 3 ) - $£(a*)) + r f x 3 (<Df (x 2 ) - $J(x 2 )) (<J>£(x 3 ) - ^(x 3 )) 



+rjx 2 ($f(x 2 ) + cDj(x 2 )) (<D£(x 3 ) + cD£ 



£ ; ei(*')^n(xi,x 3 ,x 2 ,6i,6 3 ) I .(59) 



In the annihilation diagrams, if /o is the upper meson (in the heavy b quark side), the factor- 
ization formulae are: 



"1 roo 

F a L ( a i) = 8TrC F m%f B I dx 2 dx 3 \ b 2 db 2 b 3 db 3 











x|[(x 3 - l)$£(x 3 )$/(x 2 ) - 2r^r/(x 3 - 2)$£(x 3 )$f (x 2 ) + 2r K r f x 3 <$> T K {x 3 )<5> s f {x 2 )} 
xE ai (t)h a (x 2 , 1 - x 3 , b 2 , 63) 

+ k 2 $£(x 3 )$ / (x 2 ) - 2r^r / ^(x 3 )((x 2 + l)$f(x 2 ) + (x 2 - 1)$J)] 
x J E ai (t')/ia(l-x 3 ,x 2 ,6 3 ,6 2 )l, (60) 



if P (ai) = 



-VoixC F m%f B dx 2 dx 3 b 2 db 2 b 3 db 3 lE ai (t)h a (x 2 ,l - x 3 ,b 2 ,b 3 ) x 
r K (x 3 - l)$/(x 2 ) ($&(x 3 ) + $£(^3)) + 2r/$^(x 3 )$f (x 2 )] - h a (l - x 3 , x 2 , 6 3 , 6 2 ) 
2r i ^£(x 3 )$ / (x 2 ) + r / x 2 $i(x 3 ) ($J(x 2 ) - <f>j (x 2 ))] £7 ai (0} , (61) 



LL, 



a,; = 



32irCFm B j\j2N c \ dx\dx 2 dx 3 \ b\db\b 2 db 2 § B (x\,b\) 
Jo Jo 

x^E' ai (t)h na (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) [-x 2 $£(x 3 )$/(x 2 ) 

x 2 + x 3 - l)$£(x 3 ) + (-x 2 + ^3 + 1)^N)) 



+rKrf$j(x 2 ) ((X2 - ^3 + 3)$£(x 3 ) - (x 2 + x 3 - 1)$ 

-- B ai(*'Xa(^1^2,X 3 ,6l,6 2 ) X [(x 3 - l)$£(x 3 )$/(x 2 ) 

+r^r / $£(x 3 ) ((x 2 - x 3 + l)$f (x 2 ) - (x 2 + x 3 - l)$/(x 2 )) 
+r^r/^(x 3 ) ((x 2 + x 3 - 1)$^ (x 2 ) - (1 + x 2 - x 3 )$J(x 2 



(62) 
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rl r-oo 

Ma R {ai) = 32irC F m%/y/2N c / dx x dx 2 dx 3 / b 1 db 1 b 2 db 2 ^ B (xi,b 1 ) x 

Jo Jo 

I [r K (l + X 3 )$/(X2)(^(X 3 ) - $£(x 3 )) + r /(^2 - 2)$a-(x 3 )(^(s 2 ) + (X2)) 
xE' ai (t)h na (x 1 ,x 2 , x 3 , h, b 2 ) 

- r K (x 3 - l)^ f (x 2 )(^{x 3 ) - (x 3 )) + r f x 2 <S> K (x 3 )(<$> s f (x 2 ) + $/(x 2 )) 

61,62)} • 

If /o is the lower one, 

rl roo f r 

8irC F m B f B dx 2 dx 3 6 2 d6 2 6 3 eZ6 3 | \(x 2 - l)$^(x 3 )$/(x 2 ) 

+2r^r/(x 2 - 2)*^(x 3 )*^(x 2 ) - 2r^r/x 2 $£(x 3 )<l>J(:E2) 

xE ai (t)h a (x 3 , 1 - x 2 ,6 3 ,6 2 ) + E ai (t')h a (l ~ x 2 ,x 3 ,b 2 , b 3 ) x 
x 3 <S>£(x 3 )<S> f (x 2 ) + 2r K r f $ S f (x 2 ) ((x 3 + l)$£(x 3 ) + (a* - l)*£(x 3 )) 



?LL> 



16irC F m 4 B f B dx 2 dx 3 b 2 db 2 b 3 db 3 <.E ai (t)h a (x 3 ,l - x 2 ,b 2 ,b 3 ) 
x [r f (x 2 - l)$£(x 3 ) ($/(z 2 ) + $/(x 2 )) - 2rx$^(x 3 )$/(x 2 ) 
- [2r/$i(x 3 )^(x 2 ) + r^S/O^) ($&(x 3 ) - <S>k(x 3 )) 
xE ai (t')h a (l - x 2 ,x 3 , b 2 , 63)}, 

rl roo 

— 32nC F m B / \/2~Nc I dx\dx 2 dx 3 I bidbib 3 db 3 <& B (x\,bi) 
Jo " Jo 

x ^E' ai (t)h na (xi,x 3 ,x 2 , 61, 63) x 3 <5>k{x 3 )<$> f (x 2 ) 
+r K r f ^(x 2 ) ((x 2 - x 3 + l)*£(x 3 ) - (x 2 + x 3 - l)$£(x 3 )) 
+r^r/$f (x 2 ) ((-x 2 + x 3 + 3)$£(x 3 ) + (x 2 + x 3 - l)$£(x 3 )) 

+£^(00^1^3, ^2, 61, 63) [(x 2 ~ l)^(x 3 )*/(x 2 ) 

+r^r/$J(x 2 ) ((-x 2 + x 3 + l)^(x 3 ) - (x 2 + x 3 - l)$£(x 3 ) 
+r K r / $f (x 2 ) ((^2 - x 3 - l)$£(x 3 ) + (x 2 + x 3 - l)$£(x 3 )) 



(63) 



(64) 



(65) 



(66) 



rl roo r 

M^ w {a,i) = 32vrC , F m|j/v / 2iVc J d Xl dx 2 dx 3 brfh b 3 db 3 $b(x u h)lE' ai (t) 

r f (x 2 + l)*&(x 3 )(*?(x 2 ) - $/(x 2 )) + HrO"* - 2)$ / (x 2 )($£(x 3 ) + $£(x 3 )) 
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xh na (x 1 ,x 3 ,x 2 ,b 1 ,b 3 ) - E' ai (t')h' na (xi,X3,X2,bi,b 3 ) x 
r f (x 2 - l)$£(x 3 )(<!> s f (x 3 ) - $J(x 3 )) + r K x 3 $> f {x 2 )($> p K (x 3 ) + <5> T K {x 3 ))] } . (67) 

In the above formulae, rj = mf/niB and tk = mP K /mB, tuq is the chiral scale parameter. The 
function E are defined as: 

E ei (t) = a s (t) ai (t)exp[-S B (t) - S 3 (t)}, (68) 

E' ei {t) = a 8 (t)ai(t)exp[-S B (t) - S 2 (t) - S 3 (t)]\ h=bl , (69) 

E ai (t) = a s (t) ai (t)exp[-S 2 (t) - S 3 (t)], (70) 

E' ai {t) = a s (t)a i (t)exp[-S B (t)-S 2 (t)-S 3 (t)]\ b3=b2 , (71) 

where a s is the strong coupling constant and ai is the corresponding Wilson coefficient, S is the 
Sudakov form factor. In our numerical analysis, we use the one-loop expression for the strong 
coupling constant; we use c = 0.3 for the parameter in the jet function. The explicit form of h and 
5 have been given in |18j . 
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